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Abstract 

The phase behavior of Ising spin fluids is studied in the presence of an ex- 
ternal magnetic field with the integral equation method. The calculations are 
performed on the basis of a soft mean spherical approximation using an effi- 
cient algorithm for solving the coupled set of the Ornstein-Zernike equations, 
the closure relations, and the external field constraint. The phase diagrams 
are obtained in the whole thermodynamic space including the magnetic field 
H for a wide class of Ising fluid models with various ratios R for the strengths 
of magnetic to nonmagnetic Yukawa-like interactions. The influence of vary- 
ing the inverse screening lengths Z\ and Z2, corresponding to the magnetic and 
nonmagnetic Yukawa parts of the potential, is investigated too. It is shown 
that changes in R as well as in Z\ and z% can lead to different topologies of the 
phase diagrams. In particular, depending on the value of R, the critical tem- 
perature of the liquid-gas transition either decreases monotonically, behaves 
nonmonotonically, or increases monotonically with increasing H. The para- 
ferro magnetic transition is also affected by changes in R and the screening 
lengths. At H = 0, the Ising fluid maps onto a simple model of a symmetric 
nonmagnetic binary mixture. For H — * 00, it reduces to a pure nonmagnetic 
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fluid. The results are compared with available simulations and the predictions 
of other theoretical methods. It is demonstrated, that the mean spherical ap- 
proximation appears to be more accurate compared with mean field theory, 
especially for systems with short ranged attraction potentials (when z\ and 
Z2 are large). In the Kac limit z\,Z2 — ► +0, both approaches tend to nearly 
the same results. 

PACS number(s): 64.70.Fx, 05.70.Fh, 64.60.-i, 75.50.Mm 

I. INTRODUCTION 

During the past decades, much attention has been paid to the phase behavior of ferro- 
magnetic fluid models with coupled spin and spatial interactions [1-20]. The investigations 
have been carried out using the mean field (MF) theory [1-7], the method of integral equa- 
tions [8-13], as well as Monte Carlo (MC) simulation techniques [7,8,11,14-20]. They dealt 
mainly with simplified models belonging to a so-called ideal class of spin fluids, where the 
attractive part of nonmagnetic interactions is absent (R = oo). Moreover, these models were 
considered, as a rule, in the absence of an external magnetic field (H = 0). At finite values 
of R, it has been established [1,2,4,5] that, depending on the relative strength of magnetic 
interactions, the gas, liquid, paramagnetic, and ferromagnetic states in the system may form 
phase diagrams of different topologies. For instance, an order- disorder liquid-liquid phase 
transition may appear additionally to the gas-liquid one. 

A complete picture on the phase diagram topology can only be obtained if the field is 
included (H ^ 0) in the consideration. Then the phase diagram shows the two other critical 
lines (so called wings) meeting the magnetic transition line in the H = plane at the 
tricritical liquid-liquid transition [4,5]. The gas-liquid critical point extends to a critical line 
in the magnetic field. Whether the gas-liquid or liquid-liquid critical line ends in a critical 
end point and the corresponding other critical line tends to infinite magnetic field depends 
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on the model parameters mentioned. This enlarges the number of different "phases" in the 
global phase diagram (the regions where the different topologies exist in the space of the 
microscopic model parameters, like interaction strength or range of the potentials). 

The global phase diagrams in one-dimensional i?-domain were obtained for Ising and 
Heisenberg spin fluids within MF theory [4,5]. To our knowledge, no such diagrams have 
been investigated up to now using the integral equation method. As is well recognized, 
the latter method leads to more accurate predictions. It takes into account pair correlations 
between particles in spin space, which are completely ignored by the MF approach. Previous 
integral equation studies on magnetic fluids have been restricted exclusively to ideal systems 
with Heisenberg spin interactions in the absence of an external field [8,9] or only included 
a few nonzero field values [10-13]. No integral equation calculations have been performed 
for nonideal spin fluids at H ^ 0, even within the well-known Ising model. Note that here 
we are dealing with genuine fluid models, meaning that the spatial positions of spins are 
distributed continuously, contrary to simplified lattice gas schemes [21-23], where the spins 
are positioned on fixed lattice sites. 

Due to the discrete character of spin reorientations in the Ising fluid, it can be mapped 
onto a binary nonmagnetic mixture with symmetric interparticle interactions. In this con- 
text, it should be pointed out that in recent years a lot of papers have been devoted to study 
phase properties in symmetric mixtures by means of MC simulations, the MF theory, as well 
as the Ornstein-Zernike (OZ) integral equation method [24-31]. Various closure relations, 
including the standard mean spherical approximation (MSA) [32] and a self-consistent OZ 
approach (SCOZA) [29,33], have been exploited within the latter method. Note that these 
studies considered in fact only the case when the chemical potentials of different species 
are fixed to be equal. In the language of Ising fluids this corresponds to the absence of an 
external magnetic field. 

It has been realized that despite its relative simplicity, the MSA is able to give reliable 
results for the coexistence phase boundaries including the location of critical points. The 
accuracy of the MSA scheme gets worse only when calculating critical exponents. On the 
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other hand, the more cumbersome SCOZA technique can provide us with highly precise 
results for the phase boundaries and it remains accurate even near criticality [33]. The other 
concept is the hierarchical reference theory (HRT) [30,34-36], that combines features of the 
renormalization group theory (RGT) and theoretical liquid-state approaches and allows to 
reproduce some critical exponents more precisely with respect to the MSA. For instance, 
the critical exponent f3 (which gives the curvature of the coexistence curve near the critical 
point) takes the values 7/20 and 0.345 within the SCOZA and HRT, respectively. They are 
close to the experimental and RGT prediction (3 = 0.327, contrary to the MF and MSA 
value (5—1/2 (see section III.B). 

However, the SCOZA has so far been implemented only for a restricted class of hard- 
sphere- Yukawa potentials. For these potentials, some solutions within the MSA ansatz 
can be derived [37-39] in a semianalytical form as a set of nonlinear algebraic equations 
(which should further be solved numerically). The SCOZA employs such solutions at an 
intermediate stage of the calculations. In the presence of potentials of any other structure, 
for instance, of Lennard- Jones (LJ) type or potentials with a soft-core repulsion part, the 
mathematical structure is less amenable. In addition, the SCOZA enforces the consistency 
between different thermodynamic routes in a somewhat phenomenological manner by intro- 
ducing an artificial "temperature" depending on density and concentration. The high level 
of sophistication of the SCOZA and HRT concepts leads, in turn, to substantial computa- 
tional problems, when applying them to more realistic interaction potentials. Therefore, for 
our problem we stick to a variant of the MSA, being well aware of the properties of this 
approximation near the critical point. 

In the present paper, the global phase diagram of the Ising fluid are investigated on the 
basis of the OZ integral equation method with a soft MSA closure. This allows us to obtain 
the complete thermal phase diagrams covering the whole range of the relative strength of 
magnetic interactions and other parameters of the interaction potentials. The dependencies 
of the critical temperatures and densities on the external field are analyzed in details as well. 
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II. BACKGROUND 



A. The Ising model 



The full potential energy of the Ising fluid can be written in its most general form as 



Y n N 

Z i±j i=l 



(1) 



where rj is the (three-dimensional) spatial coordinate of the i-th particle carrying spin Sj = 



±1 with i = 1, . . . , N and N being the total number of particles, r» 



rJ denotes 



the interspin separation, and H relates to the homogeneous external magnetic field. The 
exchange integral (J > 0) of ferromagnetic interactions and the attraction part (/ > 0) of 
nonmagnetic ones can be chosen in the form of Yukawa functions, 

2(z 1 a) 2 e J a 

exp[-zi(r - cr)J , 

(2) 



J(r) = 



Zi<j + 1 r 



2(z 2 af e/o- 

J ( r ) = —r exp[-z 2 {r - a)\ , 

z 2 cr + 1 r 

where Z\ and z 2 are the inverse screening lengths of the potentials, ej and ej denote the 
interaction intensities, and a relates to the particle size. The repulsion between particles 
can be modeled by a (more realistic) LJ-like soft-core (SC) potential [7], 
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rather than by the hard-sphere (HS) function 
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The multipliers 2(zi j2 (t) 2 / (zi^u + 1), entering in Eq. (2), have been used for the sake of 
convenience of comparison of our results with previous predictions (see section III). Then, 
for instance, the integrals f£° J(r)dr = 87re/<7 3 and f£° J(r)dr = 8ne ja 3 are independent of 



z\ and Z2, respectively. Within the standard hard-sphere MF theory (HSMF), such integrals 
describe the contribution of the interactions to the free energy [3]. Thus, we can say in 
advance that the HSMF results will not dependent on z\^. At z\pp = 1 (the case which 
is usually considered in theory and simulation), the multipliers go to unity, and we come 
to the usual form for I(r) and J(r). Within the soft-core MF theory (SCMF), introduced 
recently in Ref . [7] , a slight zi^-dependency should appear. Then the integrals transform to 
J °° exp(— (3(fsc){I, J}(r)dr = 87 (T, zi^nei^a 3 , where [3~ l = k^T denotes the temperature 
with k-Q being Boltzmann's constant, and j(T, z\p) is the function which takes into account 
the softness of repulsion potentials (see Eqs. (8)-(10) of Ref. [7]). For the integral equation 
approach we expect a more pronounced dependence of the results on z it2 - 

B. Integral equation approach 

1. Mapping to a symmetric binary mixture 

Since the spins Sj in an Ising fluid take only two values, ±1, we can map the system with 
N particles carrying spin 1 or —1 onto a binary mixture with N a and N b particles of type a 
and b, respectively, where N a + Nb = N. Then Eq. (1) transforms to the equivalent form 

I N a i N b N a ,N b 

U= 2 E&"(ry) + 2 EM»"(j)+ E Mrij) - HM (5) 
where M = YliLi Si = N a — N b relates to the total magnetization of the system and 

0aa(r) = Mr) = <f{r) - [I(r) + J(r)] , 

(6) 

<Pab(r) = (f) ba (r) = ip(r) - [I(r) - J{r)} 

describe the interactions between like and unlike particles in the mixture. 

In a further step, we have to rewrite the energy, the magnetic field, and the magnetization 
per particle m = M/N in terms of appropriate variables suitable for the mixture. Apart from 
the total number density p = N/V, where V denotes the volume, such variables include the 
particle concentration x and the chemical potentials. One has the concentration relations 
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In addition, we emphasize that due to a finite value of the external field term on the right- 
hand side of Eq. (5), we have to deal with a mixture prepared in an unusual way. Indeed, 
when transferring one particle from species a to species b - in the Ising liquid this amounts 
to flipping a spin from up to down and thus AM = — 2 - without changing the spatial 
coordinates, the total change in energy is equal to AU = 2H (the other terms do not 
contribute to AU due to the symmetry (fi aa = (fi bb of the particle interactions). On the other 
hand, the change in energy of the mixture is given by 

AU = AN a p a + AN b p b = p b (p, x, T) - p a (p, x, T) , (8) 

where p a and p b are the chemical potentials of species a and b, respectively. This leads to 
the identification 

Ap = fi b (p, x : T) - fi a (p, x : T) = 2H . (9) 

Condition (9) can be considered as an additional constraint imposed on the concentration 
at given values of p, T, and H, namely, x = x(p,T,H). The reason for this procedure is 
the following. In the case of a magnetic fluid, the external field (not the magnetization) is 
accessible to experiment, whereas for a mixture it is the concentration (note the chemical 
potentials). In order to study the Ising fluid in the notation of binary mixture one, therefore, 
has to fix the difference of the chemical potentials. 

2. Formulation of the integral equations 
For mixtures, the OZ integral equations have [32] the form 

Kp{r) = c aP {r) + Pi j c «7(l r - r 'l) h ipi. r> ) dr ' » ( 10 ) 

7=a,fe 

where the total h a p and direct c a p correlation functions for a pair of particles of species 
a and j3 will depend only on their separation distance, p 7 = N 1 /V is the particle number 



density of the 7-th species, and the indices a, f3 and 7 accept two values, a (spin up) and b 
(spin down). 

Eq. (10) must be complemented by an (approximate) closure relation to be solved with 
respect to h a p and c a p. The standard MSA scheme [32], proposed originally [40] for systems 
with HS repulsion (4), should be replaced in our case by the soft MSA ansatz (SMSA) 
[41,42], appropriate for a SC potential (3). It reads 

9ap(r) = exp[-/30 a/3 (r) + h a(3 (r) - c a p(r) + B a/3 (r)] , (11) 

where g a p(j') = h a p{r) + 1 denotes the radial distribution function, and 

B aP {r) = ln[l + s aP {r)\ - s aP {r) (12) 

is the bridge function with 

Sa/3(r) = h aP (r) - c af3 {r) - /3<^(r) (13) 

(no confusion may arise between the index and the Boltzmann factor (3). Formally setting 
B = in Eq. (11) leads to the hypernetted-chain (HNC) approximation [32]. 

The SMSA additionally requires the separation of the total potential a/3 in its short- 
and long-ranged parts (ff a/3 and (fi l a/3 with a/ g = (p s a/3 + <p l a/3 . There is no general procedure 
to perform such a separation uniquely for arbitrary potentials. Since the SMSA itself is 
not exact, the Yukawa potentials in the region of core repulsion r < a allow splitting to 
some extent in various ways, leading to various versions of the SMSA. Usually the splitting 
is carried out by introducing a switch function. One choice among others is to extract 
the long-ranged part using the Boltzmann-like switching exponent built on the soft-core 
potential, i.e., 

€p{r) = -[/(r) ± J(r)] exp[-^(r)] . (14) 

Such an extraction is quite natural, because for r > o the function <ft a a{r) rapidly tends 
with increasing r to the Yukawa potential —\I{r) ± J{r)] (we note that exp[— P(fi(r)] = 1 for 
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r > v^2<T, whereas lim r ^ exp[— f3(p(r)] = 0). By the replacement <p(r) — > y?Hs( r )j we come 
to the standard HS MSA with h a/3 (r) = for r < cr and c a/3 (r) = /3[J(r) ± J(r)] for r > cr. 

Another trick lies in a modification of the bridge function to the form of Eq. (12) if 
s a p(r) > and to B a/3 (r) = when s a p(r) < 0, that combines the SMSA with the HNC 
approximation. This is in the spirit of the KH closure proposed by Kovalenko and Hirata 
[43,44]. Note that the pure SMSA sometimes leads to unphysical r-domains with nega- 
tive values of g(r). The modified SMSA preserves by construction the positiveness of g(r) 
everywhere in density and temperature space. 

3. Calculation of thermodynamic quantities 

Eqs. (10) and (11) constitute a system of six nonlinear integro-algebraic equations for the 

same number of unknowns {h, c} aa , {h, c}bb, and {h, c} a b = {h, c}b a - Once the solutions are 

found, the thermodynamic quantities are calculated in a straightforward way. In particular, 

the pressure can be calculated from the virial equation of state 

9-7T a ' b r°° r\rh a 

P(p, x, T) = pk B T - ^ £ P*P? J ^9 aP (r)r 3 dr , (15) 

where p = p a + Pb- Although the energy and compressibility routes can also be used, we will 
prefer the virial route (15) because it is most easily implemented. 
The chemical potentials can be written in the form 

Pa = p* a + k B T(\ng a + 3\nA a ) (16) 

where a = a,b and A a being the de Broglie thermal wavelength (which is independent of 
p and x). Explicit expressions for the excess part of p a can be derived using the (exact) 
Kirkwood formula [45] 

p=a,b 

Here, the integration over A corresponds to the computation of the work of transferring a 
separate particle from a vacuum (A = with <p a p(r,\) = 0) to the system (A = 1 and 
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<A*/3( r ) ^) = Vafiij'))- Performing the A-integration in a manner similar to that proposed in 
Refs. [43,44], one obtains, taking into account Eqs. (10) and (11) that 



with / a S(s')ds' being equal to (1 + s) ln(l + s) - s(s + 2)/2 at s > or for s < 0. 

It is worth mentioning that since the correlation functions are calculated approximately, 
the above expressions (15) and (18) for the pressure and chemical potentials will not be ther- 
modynamically self-consistent. In particular, the virial, energy and compressibility routes 
will lead to different results. Although the differences £1X6, clS 8b rule, relatively small, they 
may distort the phase coexistence properties near a critical point. A derivation of the self- 
consistent expressions in the case of the nonideal Ising model at the presence of the external 
field is a non-trivial problem that needs a separate serious investigation. It could be solved, 
for instance, within the SCOZA by extending its present implementation from hard-sphere 
to soft-core repulsion potentials. In this paper for the sake of simplicity we will use the virial 
pressure complemented by an appropriately chosen Maxwell construction (see subsection II. 



For a general binary mixture, the densities p ' and concentrations x ' of coexisting 
phases I and II are determined at a given temperature T from the well-known mechanical 
and chemical equilibrium conditions 




(18) 



+B a/3 (r)g al 3(r) 




B(s')ds' 4vrr 2 dr , 



C). 



C. Phase separations 



P(p\x\T) 



P(p l \x u ,T), 



p a (p ,x ,T) 



p a {p ,x ,T)=p] 



i,ii 



(19) 



•a 



p b (p\x\T) 



p h {p l \x l \T) = pf l . 



In our case, they should be complemented by the condition 
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pi u - tf l = 2H (20) 

following from the external field constraint (Eq. (9)) for each phase. It is convenient to 
replace p a and //& by the sum 

H = (21) 

and the difference Ap defined in (9). Then two of the three chemical potential conditions 
in Eqs. (19) and (20) can be rewritten in the equivalent form 

Ap(p\ x\ T) = Ap(p n , x 11 , T) = 2H. (22) 

These conditions will be satisfied automatically, provided the integro-algebraic equations 
(10) and (11) are solved in conjunction with Eq. (9). Then one finds a consistent set of 
correlation functions together with the solution 

x = x(p,T,H) (23) 

for the concentration in the mixture. For the Ising fluid, the magnetization m(p,T, H) can 
easily be reproduced from x, whenever it is necessary, using relation (7). Solution (23) can 
now be inserted into the remaining conditions yielding 

P(p\T,H) = P(p u ,T,H), 

(24) 

p(p\T,H) = p(p n ,T,H), 

where the mapping from (p, x, T)-space to the new set (p, x(p, T, H),T) = (p, T, H) has been 
performed. 

Relations (24) look now like the coexistence conditions for a one-component fluid. Indeed, 
the Gibbs free energy of the mixture G = p a N a + PbNb can be rewritten in terms of iV = 
iV a + N b and M = N a - N b as G = - &h^M = pN - HM, so that the introduced 

quantity p (Eq. (21)) has the meaning of the chemical potential of the Ising fluid. It 
can be calculated using expressions (16) and (18) [which were already used when solving 
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the external field constraint (see Eqs. (9) and (23)]. Alternatively, a Maxwell-construction 
scheme has been utilized. In order to demonstrate that this scheme can be applied in its 
standard form to the magnetic fluid, let us consider the change in the Gibbs free energy. 
According to the thermodynamic relations, we have dG = —SdT + VdP + p a dN a + p b dN b = 
-SdT + VdP + pdN - HdM , where S denotes the entropy. During the isothermal (dT = 0) 
process I — > II we obtain, integrating by parts, that for the system with a fixed number of 
particles (dN = 0) in the presence of a constant external field if, the free energy increment 
is equal to AG = PV|p - J*/ 1 PdV - H(M U - Mi). On the other hand, from the definition 
of G it follows that AG = Np\Y — H(Mu — Mi). Since the increment caused by the 
change H(M U — Mi) in the magnetic energy is the same for both routes, we come to the 
Maxwell construction Q I>n = (1/p 11 - 1/p 1 ) V + jf P(p, T, H)dp/p 2 = 0, where V denotes 
the coexistence pressure. The construction guarantees that the chemical potential will be 
the same (p|{ : = 0) in both phases I and II, so that the second line of Eq. (24) transforms 
toQ(p I ,p n ,T,H) = 0. 

In such a way, the gas-liquid and liquid-liquid phase transitions of the first order can be 
determined. The second-order para-ferro magnetic transition at H = can be found as a 
boundary (Curie) curve T\{p). Below this curve, i.e., for T < T A (p), Eq. (9), being solved at 
H = with respect to the concentration x (see Eq. (23)), should have a nontrivial solution 
x = x(p, T, 0) ^ 1/2 (i.e., m ^ 0). For T > T A (p), only the trivial one should satisfy Eq. (9) 
at H = 0. Note that the trivial solution x = 1/2 (or m = 0) appears as a result of the 
symmetry of the interparticle potentials (Eq. (6)). For the same reason, the coexistence 
phases will have an identical density at concentrations x and 1 — x, or at magnetizations m 
and —m. From the structure of Eqs. (9), (15) and (18) it also follows that, if x (or m) is a 
solution to equation (9) at some value of H, then 1 — x (or —m) will automatically satisfy 
this equation at the field —H. Therefore, the phase diagrams will be symmetric with respect 
to the magnetic field H. 
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III. NUMERICAL CALCULATIONS 



A. Computational algorithm 

The set of OZ integral equations (10) was first reduced to a system of linear algebraic ones, 
h a p{k) = c a p(k) + J2 1=a ,b P~i c ot~j(k) h 1/3 (k), by applying the 3-dimensional Fourier transform 
A(k) = j v A(r) exp(ik • r)dr = / °° 4nr 2 A(r) sin(Ar)/ (Ar)dr, where A is any function of r. It 
can be presented in the compact (2x2) matrix form h(A;) = c(k)+c(k)ph(k) with p being the 
diagonal density matrix having nonzero elements [p]u — Pa = X P an d \p\22 — Pb = (1 — x)p. 

Because of the nonlinearities in the SMSA closure (Eqs. (11)-(14)) and the external field 
constraint (FC) [see Eqs. (9) and (18)], the coupled set of OZ/SMSA/FC equations have 
to be solved iteratively. The iterations have been carried out by adapting the method of 
modified direct inversion in the iterative subspace (MDIIS) [46]. At given values of p, T, 
and H the iteration starts from initial guesses for c Qj a(r) and x, and the Fourier transformed 
functions c a p(k) are calculated. Then the total correlation functions in /c-space are obtained 
analytically, h(k) = [I — c(A;)p] _1 c(/c), where I denotes the unit matrix. Applying the back- 
ward Fourier transform to h(k) yields h(r) = l/(2ir) 3 4irk 2 h(k)sm(kr)/(kr)dk. With 
the current values of c(r), h(r), and x, the residuals to the SMSA closure (11) and field 
constraint (9) are evaluated. Using them, new values of c(r) and x are updated according 
to the MDIIS corrections, and the iteration procedure is repeated until the solutions lead 
to residuals with a relative root mean square magnitude of 10~ 6 . The coexistence phase 
densities were then found by applying the Maxwell construction. 

The ratio R of the integrated strengths of magnetic to nonmagnetic interactions can be 
calculated as 

= J0 00 4ng (r)J{r)r 2 dr 

J™4irg(r)I(r)r 2 dr ' 1 ' 

where g(r) = x 2 g aa (r) + 2x(l — x)g ab (r) + (1 — x) 2 gbb(r) is the (total) radial distribution 
function of the Ising fluid. From the form of Yukawa potentials I(r) and J(r) (Eq. (2)) it 
follows that the relation (25) transforms to 
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R = - (26) 
e/ 

for z\ = z 2 , We have preferred, the latter definition even when z\ ^ z 2 , because the former 
is sensitive to the approximation made for g(r). 

The strength e appearing in the SC potential (Eq. (3)) was set to tj. This corresponds to 
a moderate softness of (p with respect to the total potential (see Eq. (1)). In the presentation 
of our results we use the dimensionless density p* = pa 3 , temperature T* = k B T/ej, external 
field H* = H/ej, and inverse screening lengths z\ 2 — z i,2&- 

B. Results for the ideal system (R = oo) 

Examples of the phase diagrams obtained within the OZ/SMSA/FC integral equation 
approach for the soft-core ideal Ising fluid with z\ — 1 at various values, H* = 0, 0.1, 0.5, 
1, 5, and oo, of the external field are shown in sets (a), (b), (c), (d), (e), and (f) of Fig. 1, 
respectively, in the (T*,p*) plane. For the purpose of comparison, the results of the SCMF 
theory and available MC simulation data [7] are also included in this figure. 

As can be seen clearly, the OZ/SMSA/FC approach leads to much more accurate predic- 
tions of the liquid-gas coexistence densities with respect to the usual version of the SCMF 
theory. Even the adjustable version, when a semiphenomenological parameter is introduced 
within the SCMF and fitted to MC data at H — > oo (see Ref. [7]), provides us with worse 
results. At the same time, the deviations between the OZ/SMSA/FC predictions and MC 
data are relatively small, especially for regions which are well below the critical point. On 
the other hand, the precision decreases when approaching the criticality point, where the 
uncertainties in critical temperature and density estimations can reach about 10-15%. More- 
over, the computations have shown that the OZ/SMSA/FC approach reproduces, like the 
MF theory, the classical values of critical exponents. For instance, the density difference in 
liquid and gas phases, p£ — p* G) appears to be proportional near a critical point to (T* — T*) /3 
independently of H* with the critical exponent (3 = 1/2 (instead of the values f3 ~ 1/3 and 
7/20 obtained within the RGT and SCOZA, respectively). 
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The influence of varying the screening length z\ of magnetic interactions on the 
OZ/SMSA/FC phase diagram is illustrated in Fig. 2. In this respect, it should be empha- 
sized that within the SCMF scheme, the results will depend on z\ very weakly (see comments 
at the end of subsection II. A). Using the more precise integral equation approach, we can 
observe an obvious z-dependence of the binodal for all values of the external field. In par- 
ticular, the dimensionless critical temperature T* and density p* increase considerably with 
rising z\ at each fixed value of H . However, the topology of the phase diagram remains 
the same and is not affected by varying the screening length. Namely, as in the case of the 
SCMF theory [7], the OZ/SMSA/FC approach predicts a tricritical point for the ideal Ising 
fluid at H* = and not a critical end point besides a gas-liquid critical point. This holds 
for z* < 5 and at least to within the relative accuracy 1CT 6 of the numerical calculations. 

The change in z\ (this quantity will be denoted below simply as z*) does not also affect 
the tendency of the critical temperature (density) to decrease (increase) monotonically with 
increasing the external field strength H*. The dependencies T*(H*) and p* c {H*) are plotted 
in detail in subsets (a) and (b) of Fig. 3, respectively. The OZ/SMSA/FC calculations 
show that in the limit of weak fields, the functions T*(H*) and p*{H*) can be cast in the 
forms T*(H*) = T c *(0) - c T (z*)(H*) 2 / 5 and p*(H*) = p*(0) + c p (z*)(H*) 2 / 5 with c T (z*) and 
c p (z*) being quantities depending only on z*. The exponent of this power law behavior is 
in accordance with the mean field tricritical exponent [4,7] (where, of course, ct and c p are 
independent of z*). With increasing H*, the functions T*(H*) and p*{H*) begin to tend 
rapidly (especially at small and moderate values of z) to their infinite field limits. For larger 
z* , the saturation regime shifts to higher values of H* . At small inverse screening length 
(z* < 0.5), the OZ/SMSA/FC and SCMF results are practically indistinguishable. 

It is worth mentioning that in the Kac limit z* — > 0, the SCMF theory should lead to 
exact results provided the equation of state of the reference system is chosen exactly too. 
Indeed, at z* — > the magnitude 2(za) 2 /(za + 1) of magnetic potential vanishes, whereas 
the screening length 1/z tends to infinity. Under such conditions, the magnetic interactions 
can be treated as an infinitesimally small perturbation to the reference potential and the 
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assumptions of the MF theory become exact. With increasing z*, the precision of the SCMF 
description goes down. Note that within the standard HSMF and SCMF theories, the 
reference system relates exclusively to nonmagnetic HS or SC repulsion. This is acceptable 
for long ranged (z* < 1) potentials. When the screening radius is short enough (z* > 2), 
a more appropriate choice of the reference system (with including a part of the Yukawa 
potential) should be performed. On the other hand, the accuracy of the OZ/SMSA/FC 
approach when evaluating T* and p* is expected to be of order of 10-15%, in a wide region 
of ^-values. This is the same accuracy as for the case z* = 1, where the direct comparison 
with the MC data could be performed (see Fig. 1). 

The OZ/SMSA/FC results for the magnetic phase transition (which takes place only at 
H* = 0) are shown in Fig. 4 (a) for the set z* = 0.5, 1, 2, 3, and 5 of the inverse screening 
length over a wide temperature and density region. It can be seen that the dependence of 
the Curie temperature T{ on p* shifts considerably to smaller values of p* with increasing 
z*. This dependence is nonlinear contrary to the HSMF prediction, where the function 
T£ = 8np* behaves linearly on p and independently of z. Within the more accurate SCMF 
and OZ/SMSA/FC approaches, the linear dependence of T£ on p* is recovered only in the 
particular (Kac) limit z* — > 0. At larger values of z* (namely, at z* > 2), the deviations from 
the limiting behavior become significant. A similar, but considerably weaker z*-dependence 
of T£ to that presented in Fig. 4 (a) for the OZ/SMSA/FC approach is observed within 
the SCMF theory, due to the existence of factor j(T*,z*) in the temperature T£ = Six^p* 
(see Ref. [7]). Note also that the difference between the SCMF and OZ/SMSA/FC functions 
T{(p*) increases with rising z*. However, even for small values of z* , where the z*-dependence 
is not so pronounced, the SCMF theory leads to worse predictions. This is demonstrated in 
Fig. 4 (b) for a particular case z* = 1 by comparison with MC results. Although here the 
both theories agree quite well, the SCMF deviations from the MC data are, nevertheless, 
slightly larger than those of the OZ/SMSA/FC theory. 
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C. Results for nonideal models (0 < R < oo) 

1. Zero magnetic field 

The phase coexistences of the nonideal Ising fluid with z[ — z\ — 1 are shown in Fig. 5 
for H* = 0, when the ratio R of strengths of the magnetic to nonmagnetic interaction is 
not too small. As can be seen, all the curves exhibit a tricritical point behavior - type I 
of the thermodynamic phase diagrams - of the same topology as the case R = oo. Note 
that because of the great number of curves, the magnetic phase transition lines have been 
omitted in Fig. 5 as well as in Figs. 6-9. They are presented in detail in Fig. 14 below. 

With further decreasing R, the shape of the phase diagrams changes in a characteristic 
way. This is illustrated in subsets (a) and (b) of Fig. 6. Beginning from the upper boundary 
value R = R u — 0.215, beside the tricritical point (TCP) a gas-liquid critical point (GLCP) 
appears in the paramagnetic phase region (at smaller densities than the tricritical point 
density) indicating that the nonmagnetic interaction is strong enough to condense here the 
system into the liquid phase. In addition, a triple point becomes visible. The TCP now 
corresponds to a liquid-liquid transition between paramagnetic and ferromagnetic phases. 
This is different from the phase behavior of type I in the region R u < R < oo, where the 
TCP relates to the transition between a paramagnetic gas and a ferromagnetic liquid phase. 
Such a new topology of the phase behavior below R u will be referred to as type II. The 
appearance of the additional critical point with decreasing R is explained by an increased 
weight of nonmagnetic attractions in the system. The nonmagnetic interaction is sufficiently 
strong to produce a gas-liquid transition before the liquid becomes ferromagnetic. 

If the nonmagnetic interaction becomes too strong, namely, when the value of R is below 
the lower boundary level Ri, i.e., R < Ri = 0.14, the TCP disappears and transforms into 
a critical end point (CEP) (see subset (b) of Fig. 6). At the same time, the GLCP remains 
and further shifts away from the CEP. Such a topology of the phase diagram in the region 
R < Ri will be defined as of type III. Note that for extremely small values of R (when 
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R <C Ri), the phase coexistence will behave like that inherent in a simple nonmagnetic fluid 
(because then I ^> J and magnetic interactions can be ignored completely). 

These three phase diagram topologies at H = have been found earlier for other systems, 
such as symmetric binary nonmagnetic mixtures [24,25,28,29], the Heisenberg fluid [4,17] or 
the Stockmayer fluid [47]. It is interesting to remark that the boundary values R u = 0.215 
and Ri = 0.14 calculated by us within the OZ/SMSA/FC approach for the Ising fluid 
correspond to 5 U = 0.646 and Si = 0.754. Here 5 = (1 — R)/(l + R) denotes the ratio of 
interparticle potentials (outside the hard or soft core) between unlike and like (see Eq. (6), 
5 is proportional to the ratio [J(r) — J(r)]/[I(r) + J(r)]) particles in the mixture. The latter 
values are very close to those (5 U = 0.65 and 5i = 0.75) reported in Ref. [29] for a symmetric 
binary mixture and evaluated within the SCOZA technique (the HRT yields S u = 0.665 [30]). 
However, the direct comparison is not possible since we used the soft core potential instead 
of the hard sphere repulsion and another value of the inverse screening length z = Z\ = z 2 . 

Varying the parameter z* can lead to a qualitative modification of the phase diagrams 
and thus to a shift of the boundaries between different topologies. This is seen in Fig. 7, 
where the cases z\ = z\ = 0.5, 1, 2, and 3 are considered at R = 0.215 (subset (a)) and 
R = 0.14 (subset (b)). From the topologies of these diagrams it can be concluded that 
the upper R u (z*) and lower Ri(z*) boundary values decrease with increasing z* = z\ = z\. 
A more complicated situation arises when z\ ^ z\ that is presented in Fig. 8. Here, the 
quantities R u and Ri should be treated as depending on both inverse screening lengths z\ and 
z\. Analyzing the set of curves in Fig. 8, it can be stated that the behavior of R u (zl, z%) and 
Ri{z*,Z2) is not monotonic. In particular, the functions R u {zl y z 2 x ) and i?^*,^) increase 
with rising z* 2 at fixed z\ = \, but they decrease with increasing z\ at constant z\ = \. 

It is worth mentioning that the above three types of the phase diagrams can also be 
observed within MF theory [4] . The disadvantage of the MF description is that it produces 
boundary values R u and Ri which are independent of z\ and z\. This corresponds, in fact, 
to the limiting behavior of R u and Ri at z\,z\ — > 0. 
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2. Nonzero magnetic field 



A set of phase diagrams for different values of the external field H and relative strength 
R of internal magnetic to nonmagnetic interactions are plotted in Fig. 9. Here, we can see 
that the change in H modifies considerably the phase coexistence curves in the nonideal 
Ising system. At large values of R, such modifications are similar to those of the ideal Ising 
fluid (compare, for example, the subset (a) of Fig. 9, R = 5, with the subset (b) of Fig. 2, 
R = oo). Since there is no magnetic phase transition at finite fields, the tricritical point 
at H* = transforms into a critical point for H* ^ (note that one has a phase diagram 
symmetric in H — > —H) and moves monotonically in the temperature-density plane with 
increasing H to the side of lower T* and higher p*. This demonstrates that the critical 
lines meet in a tricritical point. At intermediate and small values of R, the phase diagram 
modifications exhibit nonmonotonic features. For instance, at R = 0.29 (see subset (c)), the 
critical temperature, starting from the value T c (0) at H = 0, begins first to go down, reaching 
a minimum at H* ~ 2. Further it increases up to its limiting value T coo = hin^^oo T C (H), 
where T coo can be less (subset (c)) or greater (subset (d)) than T c (0). 

More complicated scenarios are observed for parameters Ri = 0.14 < R < 0.215 = R u 
(the region of topologies of type II), when the gas- liquid critical point exists simultaneously 
(from the left in p*-axis) to the tricritical one. However, only one of these two critical points 
at H = can be connected by a critical line with the critical point at H — oo. With 
increasing H, one of the two critical lines has to end. This is only possible in a critical 
end point. How this happens depends on the concrete value of R (see subsets (e), (f), 
and (g) of Fig. 9). In the infinite field limit H — > oo, the phase diagrams tend at each 
R to the gas-liquid binodals of a simple nonmagnetic fluid with the interparticle potential 
0(r) = (p(r) — I(r) — J(r). This is because then all the spins align exactly along the field 
vector, so that the product Sj • Sj will be equal to 1 (see Eq. (1)) for any pair of particles. 
A similar behavior even for finite fields can be observed for regions with too low values of 
R < Ri = 0.14 (see subset (h)), where the influence of magnetic interactions can be neglected 
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(J«J). 

The dependencies of the gas-liquid critical temperature T*^ gl \ the liquid-liquid critical 
(wing) temperature T*^ and the corresponding densities p*^ and p*^ of the nonideal 
Ising fluid on the value H* of the external field are shown in detail in Figs. 10 and 11, 
respectively. They cover the whole region of varying R and include all the three types, I 
(subset (a)), II (subsets (b) and (c)), and III (subset (d)) of the phase diagram topology. 
Note that the both phase transitions exist only in region II of the global phase diagram, 
whereas the wing line disappears for R < Ri and the gas-liquid phase transition line does not 
appear for R > R u . For type I (subset (a)), the wing line will correspond to the gas-liquid 
critical point (no liquid-liquid phase transitions are present in i?-regions corresponding to 
types I and III). As can be seen, the monotonic decrease of T*^ with rising H*, observed 
at 0.5 < R < oo, gradually transforms into a nonmonotonic function T*^ W \H*), when the 
parameter R lies in the interval [R v i,0A] with R v i = 0.196 (subset (a) of Fig. 10). The 
position of the minimum in Tf w \H*) shifts from H* ~ 3 to 1 with decreasing R. At 
R > Rvh the wing line can exist for arbitrary fields < H* < oo. For R < R vi , the wing 
line terminates at some finite value H* e (R) (subset (b) of Fig. 10) until it disappears at 
R < Ri = 0.14. 

The wing line density p*^ w \H*) also exhibits a nonmonotonic field behavior in the in- 
terval 0.6 < R < 1.25 with a maximum at H* ~ 1 to 2. Outside of this interval, it increases 
(R > 1.25) or decreases (R < 0.6) monotonically (see subsets (a) and (b) of Fig. 11). On the 
other hand, for 0.196 = R v i < R < R u = 0.215 the gas-liquid critical point ends in a critical 
end point at some finite value H* lce (R) (subset (b) of Figs. 10 and 11). For R < R vU the 
gas-liquid transition line exists for arbitrary fields < H* < oo. The critical temperature 
T* W (H*) of this transition increases always monotonically with increasing H* (subsets (c) 
and (d) of Fig. 10), whereas the critical density p^ gl \H*) is always a nonmonotonic function 
exhibiting a maximum at H* ~ 3 to H* ~ 30 depending on R (see subsets (c) and (d) of 
Fig. 11). 
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3. Van Laar point in the global phase diagram 

From the phase diagram analysis presented it follows that in the region of topologies of 
type II (Ri < R < R u ), the gas-liquid and liquid-liquid transitions can coexist simultaneously 
not only at H = but also for H ^ 0. On the other hand, it has been realized that for 
sufficiently strong fields including the limit H — > oo, we have a simple nonmagnetic-like 
phase behavior with the presence of only one gas- liquid transition (see Fig. 9). Thus, region 
II in the global phase diagram has to split into two subregions in dependence whether the 
gas-liquid (type Ha, R vt < R < R u ) or the liquid-liquid (type lib, Ri < R < R vt ) critical 
line terminates in a critical end point at some finite value of H*. The boundary in the global 
phase diagram between these two regions defines a van Laar-like point [48] . 

Our calculations have shown that this special point is identified at R v i = 0.196 with 
if t * 2.2 for p* tr ~ 0.55 and T t * ps 21.3. As can be seen in Fig. 12, there exists a some finite 
critical value H tr , where both the gas-liquid and liquid-liquid critical lines merge in one point 
in the temperature-density plane projection. It is an asymmetric tricritical point and the 
change of the critical exponent (3 from the mean field value {(3 = 1/2) to its tricritical value 
(/3 — 1/4) is seen in Fig. 12 (subset (b)). The van Laar point has been found in symmetric 
mixtures [49] but so far not in Ising liquids. In van der Waals theory the value of the van 
Laar point is R^ F = 0.279 in agreement with the value 5 = 0.564 given in [49] (note that 
A = 1 — 5 = 0.436 there). The value R v i = 0.196 is comparable to that found (at z* = 1.8) 
in [30] (their S v i = 0.67 corresponds to R v i = 0.197). In order to collect our results for the 
thermodynamic phase diagrams, the global phase diagram in i?-space is presented in Fig. 13. 

4- Magnetic critical line 

Let us consider, finally, the OZ/SMSA/FC result on the para-ferro coexistence in the 
nonideal Ising fluid at H — 0. It is shown in Fig. 14 for various values of R in a wide region of 
the temperature-density plane. In the case of the ideal fluid (R = oo), a strong dependence 
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of the Curie temperature on the screening length of the magnetic interaction (see subset 
(a) of Fig. 4) has been found. As can be observed in Fig. 14, the function T x (p) exhibits 
also an independence, especially at low densities. This is in a contrast to predictions of 
the HSMF and SCMF theories which lead to values of T£(p*) independent of R (see dashed 
curves in Fig. 14). However at larger densities (p* > 1) all the OZ/SMSA/FC functions 
begin to converge to the same curve, which is almost a straight line, and the /^-dependence 
vanishes. This line does not coincide with the HSMF and SCMF results. 

IV. CONCLUSIONS 

We have formulated a generalization of the integral equation formalism for symmetric 
binary mixtures in order to study phase coexistence properties of Ising spin fluids in the 
presence of an external magnetic field H. Mapping the spin system onto the binary mixture 
shows that the calculations for the soft-core Ising fluid at a certain magnetic field reduce in 
the mixture picture to the calculations at a certain value of the chemical potential difference 
of the constituents of the mixture. This introduces a field constraint to the OZ equations and 
modifies the MSA closure to the SMSA ansatz. It has been demonstrated that the resulting 
OZ/SMSA/FC approach is able to describe adequately the phase behavior of such models. 
Depending on the ratio R of the strengths of the magnetic and nonmagnetic interactions 
inherent in the spin fluid system, four types of thermodynamic phase diagrams have been 
identified. 

As has been established, the OZ/SMSA/FC approach provides us with more accurate 
predictions in comparison to those of the MF theory and corroborates the MF global phase 
diagram, changing of course the boundary values of R which separates the different topolo- 
gies. It is expected that some other more complicated schemes, such as the SCOZA, for 
example, should lead to a higher precision of the calculations. However, they are not yet 
developed in their present formulations to be directly applied to magnetic fluids with soft- 
core repulsion potentials. Due to their high level of sophistication, they meet considerable 
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computational difficulties in actual implementations. On the other hand, the approach pro- 
posed here can be used for systems with arbitrary potentials at relatively low computational 
costs. 

The OZ/SMSA/FC scheme can also be extended to magnetic systems where spins accept 
more than two discrete values (they map onto a nonmagnetic multi-component mixture). 
The present SMSA can be applied to a soft XY and Heisenberg fluids (replacing the dis- 
crete external field constraint (Eqs. (9) and (18)) by its continuous spin counterpart, such 
as the Lovett equation [50], for instance). These questions as well as the problem of improv- 
ing the thermodynamic self-consistency of the integral equation approach is left for future 
considerations. 
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FIGURE CAPTIONS 



FIG. 1. The liquid-gas coexistence densities p* as a function of temperature T* obtained 
within the OZ/SMSA/FC integral equation approach (bold curves) for the soft-core ideal 
(R = oo) Ising fluid with z[ — 1 and different values of the external field, H* = 0, 0.1, 0.5, 
1, 5, and oo, (subsets (a), (b), (c), (d), (e), and (f), respectively). The results of the usual 
and the adjustable versions [7] of the SCMF theory are plotted correspondingly by thin 
and dashed curves. The Gibbs ensemble MC simulation data [7] are shown as circles. The 
para-ferro magnetic phase transition (at H* = 0, subset (a)) is shown by the long dashed 
curve. 

FIG. 2. The liquid-gas coexistence densities p* as a function of temperature T* obtained 
within the OZ/SMSA/FC integral equation approach for the soft-core ideal Ising fluid at 
different values, z\ = 0.5, 1, 2, and 3, of the inverse screening length (subsets (a), (b), (c), 
and (d), respectively) as well as different values of the external field, namely, top to bottom 
and alternating solid and dashed curves, H* = 0, 0.01, 0.1, 0.3, 0.5, 1, 2, and 5 (within 
each subset). At z* = 3, an additional curve corresponding to H* — 9 is included in subset 
(d). Note that the curves corresponding to H* = oo are not shown, because they practically 
coincide with those related to H* — 5 at z* — 0.5, 1, and 2 (subsets (a), (b), and (c)) or to 
H* = 9 at z\ = 3 (subsets (d)). The para-ferro magnetic phase transition (at H* = 0) is 
plotted in the subsets by the long dashed line. 

FIG. 3. The critical temperature T* (subset (a)) and critical density p* (subset (b)) 
as functions of the external field H*, evaluated within the OZ/SMSA/FC approach for the 
ideal Ising fluid at various values of the inverse screening length z* . At z* — > 0, the result 
corresponds to the SCMF theory. 

FIG. 4. (a) The para-ferro magnetic transition temperature T£ obtained at H* = 
within the OZ/SMSA/FC integral equation approach for the soft-core ideal Ising fluid for 
different values of inverse screening length, from right to left, z* =0.5, 1, 2, 3, and 5. The 
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result of the standard HSMF theory is plotted as the dashed straight line and relates to 
the case z* -> 0. (b) The results of the SCMF and OZ/SMSA/FC theories for the case 
z* — 1 are shown as dashed and solid curves, respectively, in comparison with canonical MC 
simulation data (circles) taken from Ref. [7]. 

FIG. 5. The gas-liquid coexistence densities p* as a function of temperature T* obtained 
within the OZ/SMSA/FC approach for the soft-core nonideal Ising fluid with z\ — z\ — 1 
at H* = for larger values of R. 

FIG. 6. The coexistence densities p* as a function of temperature T* obtained within 
the OZ/SMSA/FC approach for the soft-core nonideal Ising fluid with z\ — z\ — 1 at H — 
for moderate (subset (a)) and small (subset (b)) values of R. 

FIG. 7. The coexistence densities p* as a function of temperature T* obtained within the 
OZ/SMSA/FC approach for the nonideal Ising fluid at R = 0.215 (subset (a)) and R = 0.14 
(subset (b)) corresponding to different values, z\ = z* 2 = 0.5, 1, 2, and 3, of the inverse 
screening radii. 

FIG. 8. The same as in Fig. 7 but for z\ ^ z 2 at R — 0.215 (subsets (a) and (b)) and 
R = 0.14 (subsets (c) and (d)). 

FIG. 9. The complete thermodynamic phase diagrams of the nonideal Ising fluid with 
z\ — z\ — \ evaluated using the OZ/SMSA/FC approach and projected onto the (T*,p*) 
plane at some typical values of R — 5, 1, 0.29, 0.215, 0.2, 0.19, 0.16, and 0.12 (subsets (a), 
(b), (c), (d), (e), (f), (g), and (h), respectively). The families of the diagrams in each of the 
subsets correspond to different values of the external field, H* = 0, 0.01, 0.1, 0.5, 1, 2, 3, 5, 
9, and oo (the value H* = 0.3 for R = 5 and 1, as well as H* = 20 for R = 0.12 are included 
additionally, whereas H* = 9 is excluded for R = 5 and 1). 

FIG. 10. The critical temperature of the wing lines T*( w \ subsets (a) and (b), and the 
gas-liquid critical lines T*^ 9l \ subsets (c) and (d), as a function of the external magnetic 
field H* for the nonideal Ising fluid at different values of parameter R. 
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FIG. 11. The critical density of the wing lines p*J- w \ subsets (a) and (b), and the gas- 
liquid critical lines p*J* gl \ subsets (c) and (d), as a function of the external magnetic field 
H* for the nonideal Ising fluid at different values of R. 

FIG. 12. The thermodynamic phase diagrams obtained within the OZ/SMSA/FC ap- 
proach for the soft-core nonideal Ising fluid with z\ — z\ — 1 at it! = R v i = 0.196. The phase 
coexistence curves projected onto the (T*,p*) plane are shown in subset (a) for different 
external field values, H* = 0, 0.01, 0.1, 0.5, 1, 2, 3, 5, 9, and oo. A more detailed phase 
behavior near the asymmetric tricritical point is presented in subset (b) for (bottom to top) 
H* = 1, 1.2, 1.4, 1.6, 1.8, 2, 2.2, and 2.6. Here, the gas-liquid (on the left) and wing line 
(on the right) critical points (circles) are connected by thin curves. The critical temperature 
rp*(gi,w) anc j cr j^j ca i density p*^ 9l ' w "> of the gas-liquid and liquid-liquid phase transitions are 
plotted as functions of H* in subset (c) and (d), respectively. They meet in the tricritical 
point (dots). 

FIG. 13. Global phase diagram of an Ising fluid. Type I contains at H = a tricritical 
point, where the magnetic transition line and two wing lines, existing for arbitrary ±H, 
meet. Type Ha contains a tricritical point (at H = 0) and the gas-liquid critical lines (for 
±if), ending in a critical end point at finite magnetic field, while the wing lines exists for 
arbitrary magnetic field. Type lib contains a tricritical point (at H = 0) and the gas-liquid 
critical lines (for rkH) extending to infinite magnetic field, whereas the wing lines end in 
a critical end point at finite magnetic field. Type III contains only gas-liquid critical lines 
extending to infinite values of magnetic field, and the magnetic transition line ends in a 
critical end point at H = 0. 

FIG. 14. The para-ferro magnetic coexistence curves evaluated at H — using the 
OZ/SMSA/FC theory for the soft-core nonideal Ising fluid with z\ — z\ — 1 at different 
values of the system parameter, namely, top to bottom, R = 0.1, 0.12, 0.14, 0.17, 0.215, 
0.29, 0.4, and oo. The results of the HSMF and SCMF approaches are plotted as the short- 
and long-dashed curves, respectively. 
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